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Abstract. We examine spin diffusion in a two-component homogeneous Fermi 
gas in the normal phase. Using a variational approach, analytical results are 
presented for the spin diffusion coefficient and the related spin relaxation time 
as a function of temperature and interaction strength. For low temperatures, 
strong correlation effects are included through the Landau parameters which we 
extract from Monte Carlo results. We show that the spin diffusion coefficient 
has a minimum for a temperature somewhat below the Fermi temperature with 
a value that approaches the quantum limit ~ h/m in the unitarity regime where 
m is the particle mass. We finally derive a value for the low temperature shear 
viscosity in the normal phase from the Landau parameters. 
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1. Introduction 

Cold atomic gases provide a unique possibility to study many-body physics in a 
controlled way as one can tune properties such as the strength and sign of the 
interaction as well as the population in various internal states. Fermi gases in the 
unitarity regime characterized by a diverging scattering length |a| — >■ oo are intensely 
investigated. They represent a universal realization of a strongly interacting Fermi 
gas in the sense that their properties are independent of the details of the interaction. 
Thermodynamic as well as microscopic equilibrium quantities have been measured in 
a series of impressive experiments [TJ HI |3J 21 |5] . For high temperatures, one can 
calculate their properties in a controlled way using the virial expansion [3 [HI (SI ■ For 
low temperatures, the system is strongly interacting and poses a challenging problem. 
Quantum Monte Carlo calculations are therefore presently the most reliable method 
to extract their equilibrium properties [lOl [Til [E] ■ 

There is increasing focus on the transport properties of these gases. As opposed 
to thermodynamic quantities, strong interactions can change transport coefficients by 
orders of magnitude, making them attractive to study. The shear viscosity 77 of a 
Fermi gas has been extracted from the damping of collective modes [13] as well as 
expansion experiments [14]. One can calculate the viscosity accurately for a Fermi 
gas in the unitarity limit for high temperatures |15|. Furthermore, for T <C with 
Tc the critical temperature for superfluidity, the contribution to the viscosity from 
phonons has been calculated for superfluid ^He [TB] and for a superfluid Fermi gas in 
the unitarity regime |17j . These two results agree if the sound velocity of a Fermi gas 
is inserted in the expression for ^He showing that the microscopic origin of suprfluidity 
is irrelevant [TH]. At intermediate temperatures, it is very challenging to calculate ry 
from a microscopic theory, since many effects have to be considered including particle 
self-energies [15], pairing, and vertex corrections [TSllin]- Using methods from from 
string theory, it has been shown that the viscosity of a class of strongly interacting 
many-body systems obey the bound rj/s > h/AnkB where s is the entropy density |21j . 
and it has been conjectured that this bound holds for all fluids [22]. This has made 
the study of the viscosity of cold atomic Fermi gases in the unitarity limit directly 
relevant to other fields including quark-gluon and high energy physics [23j . 

Very recently, measurements of spin transport properties of a two-component 
Fermi gas have been reported [24] . Inspired by this, we analyze in this paper the 
spin diffusion coefficient _D of a two-component homogenous Fermi gas in the normal 
phase. Using a variational approach, we calculate D both for weak and for strong 
interactions. For high temperatures, our calculations are analogous to the ones leading 
to the accurate result for the shear viscosity |15_' . In the low temperature regime where 
the system is strongly correlated, we use Fermi liquid theory to calculate D. The 
strong coupling effects are contained in the Landau parameters which are extracted 
from Monte Carlo calculations. In this way, we expect to obtain reliable results for D 
even in the unitarity limit. It follows from our high and low T results that the spin 
diffusion coefficient exhibits a minimum for a temperature below Tp . In the unitarity 
limit, we find that the minimum is D ^ h/m. Since D in general decreases with 
increasing coupling strength, this value can be interpreted as the quantum limit for 
how small D can become in atomic gases, when interactions are as strong as quantum 
mechanics allows. In analogy with the intriguing conjecture of a global minimum 
bound for 77, it would be very interesting to compare this quantum minimum for D 
with experimental results as well as with D for other strongly interacting systems. We 
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also calculate the closely related spin relaxation time, and briefly discuss how rj can 
be extracted from the Landau parameters. 

2. Spin diflFusion in the hydrodynamic regime 

We consider a gas of fcrmions of mass m in two internal states which we denote 
spin a =t, i- In equilibrium, the densities of the two components are equal, = 
ni — n/2 = kp/Gir'^, with n the total density. We shall focus on the non-equilibrium 
situation with a spatially varying magnetization Ti^(r) — n.^(r). The corresponding 
spin current j = — is in the hydrodnamic regime given by Fick's law 

j = -DVin^ - n^) (1) 

where D is the spin diffusion coefficient. A main goal of this paper is to calculate D 
as a function of coupling strength and temperature. 

In a vacuum, the interaction between atoms with spins ^ and ^ is given by the 
cross section 

dn l+p2a2- ^ ' 

with a is the s-wave scattering length and Pr the relative momentum. There are many- 
body corrections to ^ but they are small for T ^ Tp with Tp the Fermi temperature 
of the gas, even in the unitarity limit [251 [13]. For the densities and temperatures of 
interest, we can ignore the bare interaction between atoms with equal spins. For high 
temperature, the typical scattering momenta scale as ^ y/mkBT and it follows from 
([2]) that the gas is weakly interacting irrespective of the value of a. 

For low T, the gas is strongly interacting in the unitarity regime |a| — )■ 00 and there 
is presently no quantitatively reliable microscopic theory for calculating its transport 
properties. Assuming the gas is in the normal phase for T <^ Tp, we can however use 
Fermi liquid theory. Combined with the relevant Landau parameters extracted from 
non-perturbative Monte-Carlo calculations, we can in this way derive reliable results 
for D in the strongly interacting normal phase [IHl 123 ■ Since the gas is superfluid for 
T < Tc and is predicted to be a sizable fraction of Tp in the unitarity limit [121. 
it is not obvious that one can fulfill the criterion Tc < T <ti Tp for a normal Fermi 
liquid. However, we will show that the minimum of D seems to be located for T > Tc 
where pairing is not relevant. Also, one can in fact quench the superfluid order for 
T <Tchy rotating the gas pg] . 

2.1. The Landau- Boltzmann equation 

We proceed using kinetic theory for high T and Fermi liquid theory for low T to 
describe the spin dynamics of the gas. In the hydrodynamic regime, the typical 
length scale of the dynamics is much longer than the mean free path. The non- 
equilibrium distribution function is then close to a local hydrodynamic form, i.e. 
/(^(r, p) ~ l/{exp[ep — /Xcr(r)] + 1} with Cp the quasiparticle energy and n^{r) and ^^(r) 
spatially varying chemical potentials corresponding to the magnetization n^(r)— n^(r). 
Plugging the local equilibrium distributions into the left side of the steady state 
linearized Landau-Boltzmann equation for the two spin components and taking the 
difference yields [26l [29] 



/3vp-V(/it-A^;) = - ^o^\/^o) (3) 
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with Vp ~ VpEp the velocity and /"^ the distribution function in equilibrium. The 
difference in the linearized collision integrals is 

d^P2 ,^|P-P2| dal^ ^ 
(27r)3 

with m* the effective mass. We have written Sf^{r,p) — (5/^(r, p) = /°(r, p)[l — 
/"(r, p)]5'o-(r, p) with Sfa the deviation of the distribution function away from local 
equilibrium. The differential cross section for the scattering between particles with 
opposite spins is da^^/dil with ft the solid angle between the outgoing and ingoing 
relative momenta, p^ = (p4 — P3)/2 and p^ — (p2 — p)/2. There is also a term in 
Q coming from the induced interaction between parallel spin quasiparticles which we 
have not written explicitly, since it does not contribute to the variational results for 
D given below. In the following, we take the magnetization to vary in the z-direction. 
The spin current is then given by 

d^p 
(2^ 



/ 7d^('5/t-'5AK- (5) 



2.2. A variational expression for D 

The Landau-Boltzmann equation ([S]) can be written as k = with H = {I^ — 

— /°) and K — —(3vzdz{fi^ — fJ-i). Likewise, the spin current ^ can be written 
as -{K^)kBT/dz{ti^ - Hi) with the definition (A) = / d^pf°{l ~ /^A/(27r)3. Using 
this, one can in analogy with the case of charge current [33 EU] derive a variational 
bound for the spin current given by 

Here, J7 is a trial function for the deviation function By expanding U in 

polynomials, it has been shown that using the driving term in the Landau-Boltzmann 
equation as an ansatz, i.e. J7 oc k, yields very accurate results for the viscosity 77 jl9j . 
We will therefore use the ansatz U — p^ oc n appropriate for spin diffusion in the 
following. This gives 



2\2 



m' 



Xs {PzH[pz]) 



as our variational expression for the spin diffusion coefficient with Xs — d{n^ — 
n^)/9(/i^ — fii) the spin susceptibility. Equation ^ serves as the starting point 
for our results concerning spin diffusion. It is in fact the use of the ansatz U — Pz 
combined with momentum conservation, which gives that it is only the scattering 
between opposite spins that enters in the variational expression for D. 



3. High temperature limit 

We now calculate D in the classical limit T ^ Tp. Using /° — exp[— /3(p^/2to — /i)] <C 
1 with /X the chemical potential in equilibrium, the integrals in ([t]) can performed and 
we obtain 

(WV^J for T « (weak coupling) 
(t^)"^ for T » Ta (unitarity limit) 
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Here, 



= 4^a2 / dx ^"'l \ (9) 





with ksTa = fv^ /ma? is the thermal average of the scattering cross section weighted 
with p^. We have ajj- = 47ra^ and ctJ^ = 2nfi^ /mkBT in the weak coupling and 
unitarity limits respectively. To derive ([8|, we have used Xs = n/2kBT for the 
spin susceptibility in the classical limit. The temperature dependence in ([s]) can 
be understood as follows: The spin diffusion coefficient scales as D '-^ l^iV with 
^nif = 2/ncrsc the mean free path and v a typical velocity. In the classical regime, 
Zmf ~ 1/na? for weak coupling and l-^ii ^ rnksT /h^n in the unitarity limit. Since 
V ~ \fkj^Tjrn^ we recover the T^l"^ and T-^/^ scaling in (8|. The temperature scaling 
is analogous to what is found for the shear viscosity [15 . The result ([s]) for a gas in 
the classical limit T 3> TV has also recently been derived in Ref. [M] . 

4. Low temperature limit 

For T ^ Tp, the quasiparticle scattering takes place around the Fermi surface. The 
collision integral can then be reduced to 

(p.i^b.]) = ^^(fcsr)^/,,^,, (10) 

where 

r ^ sin 0(1 -cos 0) /"'"^^dfrli,! ^ 

is the cross section averaged over the Fermi surface. Here, 9 is the angle between 
the incoming scattering momenta p and P2 and cf) is the angle between the relative 
momenta Pr and p^. The cross section at the Fermi surface can be written as 

da^^ie, 



dn 16k 



[N{O)T^H0,^)f (12) 



where is the scattering matrix for opposite spins and N(Q) = m*kp/7T'^ is the 
density of states at the Fermi level. The scattering matrix can be parametrized in 
terms of the Landau parameters as [26] 

N{0)'PHO, cf>)^\ Y.^At - 3A1 + {At + AD cos0]P,(cos0) (13) 
where Pi{x) are the Legendre polynomials and 

pi 

A] = r. (14) 

' l + i^,V(2/ + l) ^ ' 



In ( 14 ) , the symmetry of the triplet part of the scattering is taken into account by 



the COS0 factor f3T]. We now assume that it is sufficient to use the I = Landau 



parameters. Using (12|-(13) in (11) and solving the resultant trigonometric integrals, 



yields when plugged into ([7| 

24/1 l + F^ [Tf 



■K^m Cl + C|/2 - C1C3 V r 



fTp\^ ( ^ ^2 for fci?a ^ 1 (weak coupling) 



X 



Snm \ T J \ 2.6 for kpa ^ 1 (unitarity limit) 




We have defined Ci = - and C3 = + and used {pD = kBTm*kl/6TT'^ 



and = iV(0)/2(l + Fg) for T/Tp < 1 to derive (15). Also, m* = m when Ff = 



For kpla] ^ 1, one has Fq = 2kFa/Ti and Fj^* ^ —2kpa/'K which gives the weak 



couphng result in (15). The spin diffusion coefficient increases as for low T. This 
is the usual effect of Fermi blocking of the scattering which enables the quasiparticles 
to travel further thereby increasing the spin current. 



To obtain the numerical value in (15 1 for kp\a\ ^ 1, we have used values for the 



Landau parameters extracted from Monte Carlo calculations for a strongly interacting 
two-component Fermi gas: Fq ~ -0.44 and Fq ~ 2 [lOl [32] giving Aq ~ -0.79 and 
A^ ~ 0.7. Since + ~ -0.1 almost fulfills the sum rule J^ii^i + ^i) = 0> it 
seems that it is reasonable only to use the I = Landau parameters. Note that the 



interaction between opposite spins ( 13 1 is attractive for these Landau parameters: 
iV(0)7~^"^ ~ —1.4 — 0.06 cos (/i. This refiects the tendency for pairing in the gas 
which suppresses the spin susceptibility by the factor 1/(1 + Fq) ~ 1/3. Since it 
is the chemical potentials and not the magnetization — which drive the Landau- 
Boltzmann equation (|3|, the spin diffusion coefficient is increased by the same factor. 

5. Minimum of the spin diffusion coefficient 

In Figs.^lpj we summarize our results by plotting in 13 as a function of T and 1/A:i?|a|. 
In Fig. [lT(a), we plot D in units of -Dhigh/(fc_Fa)^ with -Dhigh = 9'K'^^^h/32\/2m as a 
function of T in the weak coupling regime. The low and high T curves are given by 
the weak coupling hmits of (|8| an d (jlS] ). Likewise, we plot in Fig. [l] (b) £'(T)/_Dhigh 
in the unitarity limit as given by ([8| and ([l5|. Since D(T) will interpolate between 
these two limits, these results show that D(T) will exhibit a minimum. Figure [l] 
indicates that in the unitarity limit, the minimum is for T somewhat below Tp but 
above Tc ~ 0.2rF where pairing effects have to be included. The scale of the minimum 
value is -Dmin 9'K^^^h/32^/2m ~ l.lh/m. Since D decreases with increasing coupling 
strength, this result can be regarded as the quantum limit of D when the interactions 
are as strong as quantum mechanics allows. 

In Fig. [2j we plot D as a function of l/fcj?|a| in the high temperature limit 
(a) and in the low temperature limit (b) using the units Dhigh x (T/TpY^^ and 
3h/8nm x {Tp/T^ = -Diow x {Tp/TY respectively. The weak and strong coupling 



results are again obtained from (|8| and (15), and D{l/kp\a\) will interpolate between 



these two limits. We see that D decreases with increasing coupling strength kp\a\ and 
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Figure 2. The spin diffusion coefficient as a function of l/kp\a\ in the (a) high 
T and (b) low T limits. 



that it eventually flattens out in the limit kp\a\ oo. 



6. Spin relaxation time 

From our results, we now derive an expression for the spin relaxation time which gives 
the typical time between scattering events for spin dynamics. A suitable definition 
of To can be obtained from the relaxation time approximation, i.e. by assuming 
~ II ~ ^{Sf^ — Sfi)/T]j. Using this, the linearized Landau-Boltzmann equation 
^ is easily solved for 5f^ — 6f^. Plugging into ([5| then gives the usual expression 
£) = |(1 + Fq)vpTd for T <^ Tp with vp the Fermi velocity. We define the spin 



diffusion time Tp, by comparing this to (15) which gives 



9fi fTp\^ ( ^ y, for /ci?a ^ 1 (weak coupling) 



IGnkBTp \ T J [0.9 for kpa > 1 (unitarity hmit) 

for T <^ Tp. Likewise, the relaxation time approximation yields D = kBTTp)/m for 
T ^Tp. When comparing to ([s]) this gives 

TkFW [tf) for T « Ta (weak coupling) 

( ^ j for T > (unitarity limit) 

for T ^ Tp. In Fig. ([s]) we plot the high and low T limits spin relaxation time as 
given by (16) and ([Tt]). We have used the units Thigh/(fcFa)^ for the weak coupling 
case and Thigh in the unitarity limit with Thigh = 9TT^^^h/32\/2kBTp. Again, we see 
that the spin relaxation time exhibits a minimum for T < Tp in the unitarity limit. 
However, td increases monotonically with decreasing T in the weak coupling regime. 

The spin relaxation time is useful to estimate the nature of the spin dynamics 
in a particular trapped atomic gas experiment: When uto <C 1 with co the relevant 
trapping frequency, the spin dynamics is hydrodynamic whereas it is collisionless for 

UTp, » 1. 



7. Viscosity for low T 

For completeness, we briefly outline how one can obtain the viscosity rj and the viscous 
relaxation time t^ for T /Tp <^ 1 for a normal gas from the Landau parameters. The 
variational approach for calculating r/ is explained in Ref. [30]. For T/Tp <^ 1, the 




Figure 3. The spin relaxation time as a function of T in the (a) wealt and (b) 
strong coupling limits. 



viscosity is determined by the angular average of the cross section over the Fermi 



surface as in ( 10 For the viscosity however, it is the full scattering cross section 
da^^/dfl + ^daj^ /dQ which enters. Here da^J/dQ describes the induced interaction 
between parallel spins which cannot be ignored in general, even though the bare 
interaction ([2| is only between opposite spins. One has in analogy with ( 12 1 

^^j^TOrtt(.,,r (IS) 

with 

iV(O)7^^(0,0) +AnPi(cos0)cos0. (19) 

I 

Taking only the / ~ Landau parameters, we get N{0)T^'^{d, (/>) — {Af^ + Aq) cos0 = 
— 0.12 cos This means that the induced interaction between parallel spins is 
attractive but much weaker than the interaction between different spins. Following 
steps analogous to the spin diffusion case described above, we obtain 

24 T'^ 1 f T 

=-^?i;i = 0.1 — nh. (20) 

for T <C Tp- Using the relaxation time result, rj = n^^Tr,/^ to define the viscous 
relaxation time r,,, we get 

, _ 6or| 1 h _^JTfV h 

~ ^3 r2 d + |C| ksTp ^ \T ) ksTp ■ ^^'> 



Comparing (21 1 with (16), we see that the two relaxation times are qualitatively the 
same as expected. They therefore yield the same prediction for the cross-over between 
hydrodynamic and collisionless behavior. 



8. Conclusions 



Using a variational approach, we analyzed the spin diffusion coefficient and the spin 
relaxation time for a two-component homogeneous Fermi gas in the hydrodynamic 
limit. We derived analytical results in the high and the low temperature regimes 
including strong coupling effects through Landau parameters extracted from Monte 
Carlo calculations. Our results indicate that the spin diffusion coefficient exhibits a 
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minimum for a temperature below Tp but above Tc, with a value which scales as 
h/m in the unitarity regime. It would be very interesting to compare this result 
with the value of the spin diffusion coefficient in other strongly interacting systems. 
Also, one should analyze the effects of pairing on spin diffusion in atomic gases. New 
experimental insight is particularly relevant due to the lack of a controllable method 
for calculating the minimum value of D quantitatively in the strong coupling limit. We 
finally provided expressions for the shear viscosity in terms of the Landau parameters. 
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